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ABSTRACT 

We have conducted Combined Array for Research in Millimeter-wave As- 
tronomy (CARMA) observations of LiH, in absorption, toward three quasars. 
These quasars, B0218+357, PKS1830-211 and PKS0201+113, have redshifts of 
z = 0.685 — 3.387, and shift the LiH J = 1 — transition to the 1 mm and 3mm 
wavelength bands, where atmospheric absorption is sharply reduced from that 
predominating near the rest frequency of 443 GHz. We report a 3a detection 
of LiH toward B0218+357 with a column density of 1.4 x 10 12 cm -2 and place 
an upper limit on the 6 Li/ Li ratio of < 0.16. LiH was not detected toward any 
other source. 

Subject headings: astrochemistry — ISM:molecules — radio lines:ISM 



1. Introduction 

Big bang nucleosynthesis (BBN) remains our earliest probe of the universe, and plays 
a central role in post-WMAP cosmology. Light element observations and BBN theory mea- 
sure the cosmic baryon density (p# oc Vtsh 2 ); comparis on with the independent WMAP 



meas ure of Vt^h 2 fundamentally tests the hot big bang ( jSpergel et al.l 120031 ; |L arson et al. 



201l|). The result: deuterium observations spectacularly agree with WMAP, as does helium; 



(Cvburt et al. 


2003; 


Coc et al. 


2004 



Cyburt et al.l 120081 ). This discrepancy, between the CMB-predicted primordial Li and the 



stellar abundance determinations, is the primordial lithium problem. 



The lithium problem and p ossible solutions have recently been reviewed by lFieldsl (120111 ) 
and lPospelov fc Pradlerl ( 12010( 1. A serious possibility is that nonstandard processes occurred 



during primordial nucleosynth esis, altering the lithium abundance; dark matter decays can 
do this ( Uedamzik et al.ll2006f ). A signature which would test this scenario is the production 
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of the 6 Li isotope at levels 6 Li/ 7 Li ~ 0.05-0.10 flPospelovl 120071 : ICvburt et al.ll2006l ). a ratio 
orders of magnitude larger than is found in standard BB N. Indeed, recent obs ervations of 
lithium isotopes in halo stars suggest 6 Li at these levels (lAsplund et al.l 120061 ). However, 
these results are highly controversial due to the difficulty of extracting isotope ratios from 
the thermally-broadened 6707 A Li I line. 

The lithium problem may be even worse if there are additional sources of Li which 
leave traces in the halo stars where Li is observed. A very likely source of this kind is Li 
produced by the interactions of pre-Galactic or "structure formation" cosmic rays (SFCRs) 
which should be ac celerated in cosmic shocks which accompany the infall of baryonic gas onto 
protogalaxies (e.g. iMiniati et al.ll2000l ). The resulting non-primordial 6 Li and 7 Li production 



can be significant compared to their primordial levels ( IFields fc Prodanovidl2005l ) and would 
need to be subtracted from halo star abundances to determine the true primordial levels. This 
both exacerbates the lithium problem but also opens the possibility of using Li abundances in 
different cosmic environments to probe both the early universe and cosmic-ray acceleration 
in structure formation. 

Lithium also plays a role in the primordial chemis try leading to the formation of the first 
stars . Namely, the lithium hydride (LiH) molecule ( IBougleux fc Gallil Il997t iGalli fc Palla 
19981 ) contributes to the cooling of the dust-free primordial gas. Current calculations suggest 
however that the impact of LiH is not dominant due to the small primordial Li/H abun- 
dance, but the che mistry of primordial LiH is complex and a subject of ongoing interest 
teovino et aDboilh . 



Lithium thus takes a uniquely important role in the cosmology of the early universe. 
Unfortunately, lithium is difficult to measure due to its tiny abundance: in the solar system, 
Li/H ~ 10~ 9 , while in metal-poor halo stars, Li/H ~ (1-2) xlO -10 . To date, optical spec- 
troscopy of Galactic stellar atmospheres has been the only probe of Li in the early Galaxy, 
and stubborn systematic uncertainties could account for the Li problem. Any alternative 
means of observing Li in any source would be extraordinarily useful for many reasons, not 
least of which would be the different systematics. Furthermore, to date there are no stellar 
observations of extragalactic Li, and so any extragalactic information would shed qualita- 
tively new light onto this problem. 

Since the atomic lines of Li have been so challenging to observe, we have chosen to 
search for the si mplest Li containing molecu le, lithium hydride (LiH). LiH has been searched 
for previously by lCombes fc Wiklindl ( 119981 ) toward B0218+357, where the authors report a 
tentative detection in a system that is 5 km/s offset from the systemic velocity but within 
the CO velocity profile. 
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2. Sources 

The LiH molecular isotopomers have a v=0, .7=0-1 rotational trausitiou at 443.95293 



GHz ( 7 LiH) and 453.16028 GHz ( 6 LiH) f lBellini et al.l 11994! ; iPlummer. Herbst. fc De Lucia 



1984J). This region of the mm-wave spectrum is highly opaque through the atmosphere. 
However these lines are accessible in the 1 & 3 mm bands at the Combined Array for Re- 
search in Millimeter-wave Astronomy (CARMA) for extragalactic sources in redshift ranges 
of z=0. 64-1.1 and 2.85-4.2, respectively. 

Our strategy for detecting LiH in these redshift ranges is to search for interstellar LiH ab- 
sorption toward lensed quasars with high millimeter continuum flux density. The absorption 
originates from molecular clouds in the lensing galaxy along the line of sight. These galaxies 
should preferentially exhibit high optical depth in other atomic and/or molecular transitions, 
including H and CO. There are only a few candidates known to date. The three best can- 
didate lenses which meet these criteria are B0218+357, PKS1830-211, and PKS0201+113. 
Table [T] lists the source redshift, and red-shifted rest frequencies of our search lines ( 7 LiH, 
6 LiH, and 13 CO). 



2.1. B0218+357 

B0218+357 is a gravitational lens with an Einstein ring. The source is a BL Lac object 



hot spot" f 


Patnaik et al. 


Cohen et al. 


20031). The 



point sources, A & B, have a flat spectrum, while the Einstein ring has a steep spectrum. 
They are highly variable on a timescale of a few days and each is separated into a t least 



two subcomponents separated by a few mas ( iBiggs et al.l 1200 ll ; iPatnai 



A source is brighter by a factor of ~3 in the radio and both have jets (j Wiklind fc Combes 



k et al. 1995). The 



1995 



Patnaik et al 



1995 



Biggs et al.l 2001). The lensing source is a spiral galaxy at a redshift of z~0.685 



cm H line (jCarilli et al. 



Wikl ind &: Combeslll9951). Atomic and molecular li nes including the 21 



1993h. HoCO dMenten fc Reidlll996h. HoO. CS JCombes fc Wiklind 



1997 ). CO, HCO+, HCN Jwiklind fc Combeslll995j) . NH, jHenkel et al.ll2005h. OH fcanekar et al 
2003), H 2 CO anti inversion dZeiger fc DarlineJboioh . and LiH (tentative) JCombes fc Wiklind 
19981 ) have been detected toward the A component. The line center of the LiH tentative 
detecti on is off by ~5 km/s com pared to other detected lines, but within the CO velocity 
profile dCombes fc Wiklind! Il998l ). 
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Table 1. Frequencies of Observed Lines a 

Source z 7 LiH b 6 LiH c 13 CO d Ref 

B0218+357 0.68466 263.526648 268.992129 261.634494 1 

PKS1830-211 0.88582 235.416319 240.298798 233.726001 2 

PKS0201+113 3.387144 101.194035 103.292775 100.467449 3 

a All frequency units are GHz. 

b The unshifted rest frequency of the J = — 1 7 LiH transition is 
443.95280 GHz 

c The unshifted rest frequency of the J = — 1 6 LiH transition is 
453.16028 GHz 

d The unshifted rest frequency of the J = 3 — 4 13 CO transition is 
440.7651668 GHz 



Ref erences. (1) ICombes fc Wiklindl (119981 ): (2) 
fl2008h : (3) iKanekar et all koQ'A ~ 



Muller & Guelin 
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2.2. PKS1830-211 



PKS1830-211 was one of the first quasars that was found to be lensed. The source 
blazar is at a redshift of z~2.5 and appears a s two bright images in a fain t Einstein ring at 



sub-arcsecond resolution (ILidman et al.lll999t ISubrahmanyan et al.lll990f ). The two images 



are NW and SE of the lens cente r and are separated by ~1"; the NW component is the 



(Jin et al. 


2003; 


Nair et al. 


2005) 



cause notable changes in absorption line features in the lensing galaxy (IMuller fc Guelin 



2008 ). The lens galaxy is a face on spiral at a redshift of z~0.886 (IWinn et al. 



2002 



Wiklind fc Combed 119981 ) . One of the spiral arms lies along the line of sight to the S W 
image and is the primary source of molecular and atomic absorp tion lines, although a small 
percentage also originates from the NW source (IWinn et al.ll2002l ). Nume rous molecular lines 



have been detect ed toward PKS1830-211 including CO, HCN, HCO+ ( Wik 



1995|), CS. HoO (ICombes fc Wiklindlll997h. NoH+ (Wiklind fc Combeslll998f). CoH. HC N. 



C 3 H 2 (IMenten et al.lll999h . NH 3 (IHenkel et alJ 12005( 1. and H 2 CO (IMenten fc Reidlll996h . 



ind fc Combes 



2.3. PKS0201+113 



PKS0201+11 3 is a radio loud quasar at a redshift of z~3.61 ( ICondon et al.l Il977 



White et al.lll993l ). VLBI observation s indicate that its linear size is small <40 pc and has a 
small er secondary component 2"south ( Hodges et al.|[l984 ; Morabito et al.Ml986 ; Stanghellini et al 
19901 ). The secondary source emits on 1% of the total system flux (at 22 GHz). PKS0201+113 
shows damp ed Lyman a along the line of sight, several atomic lines , and 21 cm absorption 
at z~3.387 Jwhite et al.lll993l : bva et al.lll998l : Kanekar et al.ll20Q7h . 



3. Observations 

All observations were conducted with the Combined Array for Research in Millimeter- 
wave Astronomy (CARMA). In all cases the source is unresolved in our beams, and thus 
appears as a point source. PKS0201 + 113 was observed in 2008 February when CARMA 
was in its B configuration giving a typical synthesized beam of 0.88x0.70". The u — v 
coverage of the observations gives projected baselines of 25.2-319.6 kA (74.6-946.9 m). All 
three correlator bands were set up with 62 MHz bandwidths and positioned with overlapping 
edge channels. This gives full frequency coverage across the J — 0-1 7 LiH line with a velocity 
resolution of ~2.9 km/s. Uranus was used as a flux density calibrator, the internal noise 
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source was used to correct the passbands of each window, and 0238+166 was used to calibrate 
the complex antenna-based gains. The source was then self-calibrated to refine the solution . 



The data were calibrated and imaged using the MIRIAD software package (jSault et al.lll995l ). 
We obtained about ten hours of good data on this source. 

PKS 1830-211 was observed in 2008 March and April, when CARMA was in its D config- 
uration giving a typical synthesized beam of 3.6x2.1". The u-v coverage of the observations 
gives projected baselines of 3.79-95.03 kA (4.8-121.1 m). All three correlator bands were 
set up with 62 MHz bandwidths, giving a resolution of ~1.2 km/s. Two of the bands were 
configured to overlap at the edges to cover the J = 0-l 7 LiH line in the lower sideband and J 
= 0-1 6 LiH in the upper sideband. The remaining band was configured to observe the J = 
3-4 13 CO line in the lower sideband, this choice was driven by the frequency configuration. 
Neptune and MWC349 were used as flux density calibrators, 3C273 and 3C454.3 were used 
to correct the passbands of each window, and 1733-130 and 1911-201 were used to calibrate 
the antenna based gains. The source was then self-calibrated to refine the solution. We 
obtained only 4.5 hours of good data due to poor weather. 

B02 18+357 was observed in 2008 August and September, when CARMA was also in its 
D configuration, giving a typical synthesized beam of 1.7x1.5". The u — v coverage of the 
observations gave projected baselines of 6.7-129.2 kA (7.7-147.1 m). Two of the correlator 
bands were set up in 62 MHz mode, giving ~1.1 km/s resolution, and the third was in 31 
MHz mode, giving ~0.55 km/s resolution. Band 1 in the lower sideband was centered on the 
J = 3-4 13 CO line. Bands 2 and 3 in the lower sideband were set to observe the J = 0-1 7 LiH 
transition. Band 1 in the upper sideband was set to observe the J = 0-1 6 LiH transition. 
The upper sideband of the other two bands were for continuum. Uranus was used as the 
flux density calibrator, 3C84, 3C454.3, and the internal noise source were used to correct the 
passbands of each window. The data were self-calibrated for the antenna based gains. We 
only obtained about seven hours of good data on this source due to poor weather. 



4. RESULTS & DISCUSSION 
4.1. B0218+357 

Figure [1] shows our continuum map of B0218+357. The two point sources (A and B) 
and the "hot spot" are labeled and the synthesized beam is in the lower left corner. The 
contour levels are ±12er, ±24<t, ±36cr, a = 4.5 mJy/beam and the peak flux is 885 
mJy/beam. Note that we do not detect the "hot spot" (even at our 3 a cutoff). This 
peak flux value is within the uncertainties of the expected value at our rest frequency from 
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Figure 2 of (jCombes Wiklindl 119971 ). Figure [2] shows our 3a detection of 7 LiH. The noise 
level of the spectrum is denoted by the error bar and the dashed line denotes the expected 
rest velocity of the line. This line was detected toward the A component, where all other 
line detection have been made to date and is at the expected rest velo city. For absorption 
features the total beam a veraged column density can be calculated from (ILucas fc Lisztlll993l ; 
Greaves fc Nymanl 119961 ) 



where Q r is the rotational partition function, E\ is the lower state energy, Sy? is the product 
of the line strength and dipole moment, and the opacity, r, is defined below 



MB 



[l-e- T )(J(T ex )-J(CB)-T ( 



c 



(2) 



giving 



-In 



1 - 



Ti 



MB 



(3) 



J(T ex )-J(CB)-T c \ ' 

wher e Tmb is the main beam brightness temperature of the line, J(T) is defined as flRohlfs fc Wilson 
2000h 

J(CB) is the cosmic microwave background (2.73(1 + 2;)), and Tc is the continuum brightness 
temperature. For o bservations not done i n a temperature scale one can convert from Jy 
beam" 1 to K using faohlfs fc Wilsonlhoooh 

1.22/q 



K, 



(5) 



where 1$ the is peak intensity of the line in Jy beam -1 and B is the beam filling factor 
(incorporating both the continuum and absorption source factors). In all equations the 
excitation temperature T ex is assumed to be 10 K. 

From these equations we get an opacity (r) o f 1.7 and a total column de nsity of 1.4 x 10 12 
cm -2 for 7 LiH. This value is similar to that of ( Combes fc Wiklind 1997). Calculated H 



column densities for this source range from 5 x 10 21 - 5 x 10 23 cm 2 (IMenten fc Reidlll996 



Combes fc Wiklindl Il995l ) . From these values we calculate a Li H/H? ratio of 2.9 x 10 
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2.9 x 10~ 12 . If we assume that the gas is cold and dark, as did ICombes fc Wiklindl (11998I ). 
we have f(H2) = 0.5. T he Li/H ratio for a sourc e with the redshift of B0218+357 can be 
estimated to be ~ 10~ 9 (jCombes fc Wiklind 1998 ). From these values we calculate a LiH/Li 
ratio of 0.287 - 0.003. 
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Figure [3] shows our non-detections toward B0218+357. a) shows 13 CO and b) shows 
6 LiH. Using equation [TJ we can set la upper limits of 4.8 x 10 15 cm -2 and 2.3 x 10 11 cm -2 , 
respectively. From this we c an calculate a 6 Li/ 7 Li upp er limit of <0.16. The J = 2 — 1 13 CO 
transition was detected by (jCombes Wiklindl Il995l ) . thus given our non-detection we can 
set limits on the temperature. Given our 1 a upper limit we can only set a temperature 
upper limit of 75 K. 
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Fig. 1. — B0218+357 continuum map. The contours are ±12er, ±24a, ±36<r, a = 
4.5 mJy/beam. The synthesized beam is plotted in the lower left corner and the known 
continuum sources (A and B) and the 'hot spot' (HS) are labeled. 
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Fig. 2. — LiH detection toward B0218+357. The noise level of the spectrum is denoted by 
the error bar in the lower left corner. The dotted line denotes the expected rest velocity of 
the source. 
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Fig. 3. — Non detections toward B0218+357. The dashed line denotes the expected rest 
velocity of the lines and the error bar gives the rms noise level for each spectra, a) 13 CO non 
detection b) 6 LiH non detection. 
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4.2. PKS1830-211 



Figure H] shows our continuum map of PKS1830-211. The two point sources (NE and 
SW) are labeled. The synthesized beam is in the lower left corner and the contour levels are 
±12cr, ±24<t, ±36<x, a = 9.0 mJy/beam and the peak flux is 953 mJy/beam. Figure [5] 
shows our 3.4a detection of 13 CO toward the SW component of PKS1830-211. The error 
bar denotes the rms noise level of the spectrum. The grey area indicates the velocities over 



which absorption peaks have been detected previously (e.g. iMuller et al.ll2006l ; iMenten et al. 



20081 ). The arrow indicates the line center. In order to increase the S/N of our spectra we 
inverted the data by averaging adjacent together, decreasing the velocity resolution to 2.5 
km/s. 

From equation [1] we calculate a column density of 2. 8 x 10* 6 cm -2 for 13 CO. Comparing 
this value with the 12 CO value from IWiklind fc Combed (119981 ) gives a 12 CO/ 13 CO ratio of 
~72. While this value is a bit high compared to most Galactic values it is not unreasonable, 
since 13 C is formed primarily in low and medium mass stars during their red giant phase, 
and there are fewer of these at highe r z, one could reasonably expect a higher ratio toward 
PKS1830-211 Jwilson fc Roodlll994[ and references therein). 



Figure shows our non-detections toward PKS1830-211. a) shows 7 LiH and b) shows 
I. Using 
respectively. 



6 LiH. Using equation [T] we can set la upper limits of 3.4 x 10 10 cm 2 and 3.7 x 10 10 cm 2 , 



4.3. PKS0201+113 

Figure [7| shows our continuum map of PKS0201+113. The synthesized beam is plotted 
in the lower left corner and the contour levels are ±12(7, ±24er, ±36a, a = 926 ^Jy/beam 
and the peak flux is 169 mJy/beam. 

We did not detect 7 LiH toward this source. Figure [8] shows the spectrum of the non- 
detection. The error bar denotes the la rms noise level. From equation [1] we calculate an 
upper limit to the column density of 1.0 x 10 10 cm -2 . 
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Fig. 4. — PKS1830-211 continuum map. The contours are ±12<r, ±24o", ±36cr, a = 5.6 
mJy/beam. The two point sources (NW and SE) are labeled. The synthesized beam is in 
the lower left corner. 
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Fig. 5. — 13 CO toward PKS1830-211. The error bar denotes the la nois e level of the s pectr a 
and the grey region indicates the range of expected rest velocities from iMuller et al.l ( 120061 ) ; 
Menten et all (|2008h . 
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Fig. 6. — Non detections to ward PKS1830-211. The grey regi o n ind icates the range of 
expected rest velocities from IMuller et al.l ( 120061 ) ; iMenten et al.l ( 120081 ) and the error bar 
gives the rms noise level for each spectra, a) 7 LiH non detection b) 6 LiH non detection. 
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Fig. 7.— PKS0201+113 continuum map. The contours are ±12a, ±24a, ±36a, a = 926 
/iJy/beam. The synthesized beam is plotted in the lower left corner. 
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Fig. 8. — Non detection of 7 LiH toward PKS0201+11 3. The grey region indicates the range 
of expected rest velocities from iKanekar et all (120071 ) and the error bar denotes the la rms 
noise level. 
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5. Implications and Conclusions 

We have presented CARMA observations of the J = — 1 transitions of 7 LiH and 6 LiH, 
and the J = 3 — 4 transition of 13 C0 toward three distant quasars, B0218+357, PKS1830- 
211, and PKS0201+113. While these lines were not detected in most of the sources, due 
primarily to the lack of sensitivity, some were detected. The 7 LiH line was detected toward 
only one source, B0218+357, at a ~3-a level. From these observations we calculate a beam 
averaged column density of 1.4 x 10 12 cm~ 2 . We also searched for the J = — 1 transition 
of 6 LiH toward B0218+357, but did not detect it. We were able to set an upper limit on the 
6 Li/ 7 Li ratio of <0.28 for this source. Unfortunately, this upper limit does not indicate which 
model for BBN is correct. Further observations to detect the 6 LiH transition are necessary 
to resolve this. 

Ultimately one wants to know the primordial Li abundance. Observing and detecting 
LiH is just the first step in this process. The next step is determining the Li column density 
from the LiH. Doing this requires detailed knowledge of the formation and destruction pro- 
cesses of LiH, both in the gas phase (for both high and low z objects) and on grain surfaces 
(for low z objects). To date this knowledge is notably incomplete, especially over the range 
in z in which we can observe LiH. Additionally the H column density needs to be determined 
to good accuracy. Recent determinations of the H column density in many of these sources is 
based on CO/H ratios, using CO lines that may be saturated, and assuming that the CO/H 
ratio is constant for all z, which is not likely to be the case. Thus, while LiH offers a new 
probe of extragalactic lithium, much work is still needed before LiH yields the primordial 
abundance. 
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